Microbi 




ORIGINAL RESEARCH 

Identification of a metagenomic gene cluster containing a 
new class A beta-lactamase and toxin-antitoxin systems 

Ken Vercammen 1 , Tamara Garcia-Armisen 2 , Nathalie Goeders 3 , Laurence Van Melderen 3 , 
Josselin Bodilis 4 & Pierre Cornelis 1 

1 Department of Bioengineering Sciences, Research group Microbiology and VIB Department of Structural Biology, Vrije Universiteit Brussel, 
Pleinlaan 2, 1050 Brussels, Belgium 

2 Ecologie des Systemes Aquatiques, Universite Libre de Bruxelles, Campus de la Plaine, CP 221 Boulevard du Triomphe, 1050 Brussels, Belgium 
3 Genetique et Physiologie Bacterienne, Universite Libre de Bruxelles, IBMM-DBM, 12 rue des Professeurs Jeneer et Brachet, 6041 Gosselies, 
Belgium 

4 Microbiologie Environnementale et Biologie Evolutive, UFR Sciences et Techniques, Universite de Rouen, Place Emile Blondel, 76821 
Mont-St-Aignan, France 



Keywords 

Antibiotic resistance, metagenomic DNA, 
toxin-antitoxin systems, /J-lactamase. 

Correspondence 

Pierre Cornelis, Department of 
Bioengineering Sciences, Research group 
Microbiology and VIB Department of 
Structural Biology, Vrije Universiteit Brussel, 
Pleinlaan 2, 1050 Brussels, Belgium. 
Tel: +3226291906; Fax: +3226291902; 
E-mail: pcornel@vub.ac.be 

Funding Information 

Nathalie Goeders has a FRIA fellowship. 



Received: 30 January 2013; Revised: 15 May 
2013; Accepted: 6 June 2013 

MicrobiologyOpen 2013; 2(4): 674-683 



doi: 10.1002/mbo3.104 



Abstract 

Several reports mention the presence of antibiotic resistance genes in natural 
and polluted environments, but many studies are based on their detection via 
polymerase chain reaction (PCR amplification of known genes and not on an 
activity screening. We constructed a metagenomic fosmid bank from DNA iso- 
lated from a polluted river in Brussels, Belgium, the Zenne. A total of 120,000 
clones were pooled and plated directly on solid media containing different anti- 
biotics. Several clones were isolated which could grow in the presence of ampi- 
cillin. The DNA from several clones was extracted and subjected to restriction 
analysis and, based on their restriction pattern, two different clones were found. 
One of the clones was selected for further study as it showed a higher level of 
resistance to different /^-lactams antibiotics (ticarcilline and ceftazidime). To 
find out which gene is responsible for the resistance, an in vitro transposon 
mutagenesis was performed and clones having lost the resistance phenotype 
were analyzed via inverse PCR amplification. Several clones had an insert in a 
gene encoding a new type of /^-lactamase. The amplified fosmid DNA was fully 
sequenced revealing an insert of 41 kb containing 39 open reading frames 
(ORFs). Transposon insertions inactivating the resistance to /^-lactams were also 
found in the ORF upstream of the blaA gene, encoding an aminotransferase, 
suggesting a polar effect on the transcription of the gene downstream. In addi- 
tion, other genes were found such as histidine biosynthesis genes, which were 
found to be scattered on the insert, a relA/spoT gene, and genes belonging to 
type II toxin-antitoxin system. This predicted system was experimentally 
validated in Escherichia coli using an inducible expression system. 



Introduction 

It is now well established that natural environments can 
be a reservoir of antibiotic resistance genes (Martinez 
et al. 2009; Allen et al. 2010). The increased use of antibi- 
otics in agriculture and in medicine has also caused a 
massive release of antibiotics in streams and rivers, allow- 
ing the selection of antibiotic resistant bacteria (Baquero 
et al. 2008; Martinez et al. 2009; Allen et al. 2010). As a 
result of industrial activity, toxic metals are also released 



in streams selecting for bacteria able to resist both metal 
and antibiotic imposed stresses (Calomiris et al. 1984; 
Baker- Austin et al. 2006). A recent study revealed a posi- 
tive correlation between the levels of metals and resistance 
genes (Knapp et al. 2011). For metals, as well as for anti- 
biotics, efflux via efflux pumps is involved in resistance 
(Grass et al. 2011; Kim et al. 2011). Some metals such as 
copper or chromium can cause an important production 
of reactive oxygen species (ROS) via a Fenton reaction 
(Lloyd and Phillips 1999), and recently, it has been pro- 
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posed that many antibiotics also generate an oxidative 
stress (Kohanski et al. 2010; Martinez and Rojo 2011). 
Metal and antibiotic resistance genes can be moved 
between bacteria via horizontal transfer of mobile genetic 
elements (de la Cruz and Davies 2000; Gootz 2010; 
Andam et al. 2011; Zhang et al. 2011). Two approaches 
can be used to detect antibiotic resistance genes from me- 
tagenomic DNA. The first one is based on the amplifica- 
tion of known antibiotic resistance genes using gene- 
specific markers, whereas the second, termed functional 
metagenomics, involves the construction of metagenomic 
libraries, which are further screened for expression of drug 
resistance (Riesenfeld et al. 2004; Allen et al. 2010). With 
the power of new sequencing methods, it is also possible 
to sequence the inserts of metagenomic clones, hence dis- 
covering the genomic context. In this study, we present 
the construction of a polluted stream metagenomic 
library, which was functionally screened for beta-lactam 
resistance, and the identification of a new beta-lactamase 
gene as well as its genomic context. We also present evi- 
dence for the presence and biological activity of toxin- 
antitoxin genes. 

Materials and Methods 

Strains and plasmid constructions 

See Table 1 for the list of used strains and plasmids. 

Media and growth conditions 

Yeasts were grown on 863 media (1% yeast extract, 1% 
bactopeptone, and 2% glucose). Luria-Bertani medium 



(LB) and M9 minimal medium (KH 2 P0 4 (22 mmol/L), 
Na 2 HP0 4 (42 mmol/L), NH 4 C1 (19 mmol/L), MgS0 4 
(1 mmol/L), CaCl 2 (0.1 mmol/L), NaCl (9 mmol/L), vitamin 
Bl (1 jUg/mL)) supplemented with casamino acids and car- 
bon sources were used to grow bacteria. Ampicillin and chl- 
oramphenicol were added at respective final concentration of 
50 and 12.5 /ig/mL. Isopropyl /?-D-l-thiogalactopyranoside 
(IPTG) was used at 0.1 mmol/L. 

Freshwater sampling and treatment 

Freshwater was collected in November 2008 and February 
2009 from the surface of the Zenne river at the Buda 
bridge (Haren, Belgium) situated at an industrial site. The 
Zenne water (2 L) was transferred into sterile bottles at 
4°C and immediately filtered using a 1.2- (im glass fiber 
filter (Whatman, GE Healthcare) to eliminate the sedi- 
ments and further through a 0.2-^im filter to capture the 
microorganisms. The 0.2 fim filters with the retained 
microorganisms were stored at — 80°C and thawed at 
room temperature just before extraction of the DNA. 

DNA extraction 

DNA from the filters was extracted using the Metagenom- 
ic DNA Isolation kit for Water from Epicentre®. The final 
concentration of DNA was 0.27 /.(g/mL. 

Metagenomic library construction 

Electrocompetent Escherichia coli TransformMax EPI300 
T1R cells (Epicentre®, Madison, WI) were used for library 
construction and grown on LB broth at 37°C containing 



Table 1. Strains and plasmids used for this study. 



Strain or plasmid Genotype or plasmid properties Reference or source 



Yeast 

23344c 
Bacterial strains 

MC1061 

MG16B5 

DJ624 Aara 
Plasmids 

pBAD33 

pBAD33lev 

pBADIev-relE93 

pBADIev-ORF24 

pBADIev-ORF26 

pKK223-3lev 

pKK223-3lev-relB93 

pKK223-3lev-ORF2S 
pKK223-3lev-ORF27 



MAT ura3 

araD139,D(/ac/ POZY)74, galU, galK, rpsL, A(araA BC-/eu)7697, hsdR, mcrB, Sm s 
Wild-type E. coli K-12 
MG1655 Mac malPy.laclq 

P15A, CamR, pBAD promoter 

pBAD33lev derivative containing URA3 selection gene and a yeast replication origin 
pBAD33lev derivative containing the gene relE 93 under the pBAD promoter 
pBAD33lev derivative containing the ORF24 under the pBAD promoter 
pBAD33lev derivative containing the ORF26 under the pBAD promoter 
ColE1, AmpR, pTac promoter 

pKK223-3lev derivative containing the relB g3 gene under the control of 
the pTac promoter 

pKK223-3lev derivative containing ORF25 under the control of the pTac promoter 
pKK223-3lev derivative containing ORF27 under the control of the pTac promoter 



Our laboratory 

Casadaban and Cohen (1980) 
Xiao et al. (1991) 
Our laboratory 

Guzman et al. (1995) 
This work 
This work 
This work 
This work 

Brosius and Holy (1984) 
This work 

This work 
This work 
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chloramphenicol (Cm) (12.5 /.(g/mL). The metagenomic 
library was constructed using the CopyControl Fosmid 
Library Production Kit from Epicentre®. Briefly, the 
extracted DNA was blunt ended at the 5'-P0 4 - site, 
ligated into the pCClFOS vector and packaged into X 
phages using MaxPlax™ Lambda packaging extracts (Epi- 
centre®). E. coli EPI-300 T1R cells were infected with the 
packaged bacteriophages and the clones inoculated on LB 
plates containing 12.5 /(g/mL Cm. Recombinant clones 
growing on Cm plates were scraped from the selective 
agar plates into 10 mL of LB medium plus 20% glycerol 
and stored in aliquoted 1 mL pools at — 80°C. The final 
titer was 10 n colony- forming units (CFU) per mL. 

Control of the quality of the metagenomic 
library 

Twenty colonies growing on selective medium plus Cm were 
selected randomly and grown in 1 mL LB medium plus 
Cm plus the Copy Control Autoinduction Solution from 
Epicentre. The DNA was extracted according to the manu- 
facturer protocol (Epicentre) and was restricted using the 
Fast Digest Restriction enzymes Ndel, Pstl, and Hindlll 
(Fermentas, Thermo Scientific). 

Functional screening of the metagenomic 
library 

The E. coli TransformMax EPI300 TlR::pCClFOS cells were 
tested for ampicillin susceptibility using the minimum 
inhibitory concentration assay (MIC, Oxoid, Thermo Scien- 
tific) and were found to be sensitive to less than 10 /(g mL 
Amp. The metagenomic library was therefore plated on LB 
agar containing 12.5 /(g/mL Cm and 50 /jg/mL ampicillin. 
The ampicillin resistant clones (20 clones) were also tested 
for resistance to other /Mactams, including amoxicillin, 
aztreonam, cefepime, ceftazidime, ceftriaxone, meropenem, 
piperacillin, and ticarcillin, using the disk assay (each disk 
contains 30 /ig of one antibiotic). The resistant clones were 
analyzed by restriction digestion using fast digest Pstl and 
Ndel restriction enzymes (Thermo scientific). 

In vitro transposon mutagenesis 

The EZ-Tn5<KAN-2> insertion kit (Epicentre®) was used for 
the generation of ampicillin-sensitive clones. Briefly, the fosmid 
carrying the ampicillin resistance gene was purified using the 
Qiagen miniprep kit. First a 1:1000 diluted overnight culture 
was induced for 5 h at 37°C at 225 rpm with the CopyControl 
Autoinduction Solution (Epicentre®) in order to increase the 
plasmid copy number. The purified recombinant fosmid was 
then incubated for 2 h with the transposon carrying the 
kanamycin resistance gene and the transposase under the 



conditions recommended by the manufacturer. This mixture 
was used for transformation of the E. coli TransformMax 
EPI300 T1R cells. The transposon mutants were grown on LB 
containing Cm and LB containing Cm and Amp (50 /.(g/mL). 
The fosmids of the candidates having lost their resistance to 
ampicillin were extracted, as mentioned above, and sent for 
sequencing using the KAN-2 FP-1 forward primer (5' AC- 
CTACAACAAAGCTCTCATCAACC 3') and KAN-2 RP-1 
reverse primer (5' GCAATGTAACATCAGAGATTTTGAG 3'). 

454 pyrosequencing and gene 
annotation 

The clone containing the multiresistant ampicillin fosmid 
was sent for sequencing. Briefly, the Amp R clone was grown 
overnight at 37°C and purified as mentioned above followed 
by an 1-hour incubation with Riboshredder RnaseA (Epi- 
centre). The 454 pyrosequencing reaction was done by 
DNAvision (Gosselies, Belgium). The sequences were assem- 
bled using CAP3 version 3 (available at Mobyle Pasteur: 
http://mobyle.pasteur.fr/cgi-bin/portal.py#forms::cap3) and 
annotated by the free online program Softberry (http://li- 
nuxl.softberry.com/berry.phtml?topic=fgenesb&group=pro- 
grams&subgro up =gfindb ) . 

Phylogenetic analysis 

In order to get a broad picture of the diversity of the 
Class A Beta-lactamase, 29 reference sequences were cho- 
sen from the study of Hall and Barlow (2004). In addi- 
tion, the five best BLAST hits were added with the new 
sequences identified in this study. Altogether, the 35 pro- 
tein sequences were aligned with ClustalW2 with default 
parameters (Larkin et al. 2007) implemented in Seaview 
4.2.4 (Gouy et al. 2010). After removing positions with 
gaps, conserved blocks of the multiple alignments were 
selected by using Gblocks 0.91b with default parameters 
(Talavera and Castresana 2007). A neighbor-joining tree 
was inferred using MEGA v5.0 (Tamura et al. 2011) with 
Poisson correction. The phylogeny was rooted with the 
group that includes CGA1, CME2, and CFXA2 because 
those genes are in the chromosomes of species within the 
Cytophaga-Flexibacter-Bacteroides (CFB) group, the group 
that is closest to the root of the Eubacterial tree among 
the groups represented in this phylogeny (Hall and Bar- 
low 2004). The node supports (on the branches) were 
evaluated based on 1000 bootstrap replicates and the bar 
indicates sequence divergence. 

Cloning of toxin and antitoxin genes 

Primers used to amplify genes are flanked by 40 bp 
homologous to the cloning sites of the pBAD33-lev or the 
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pKK223-3-lev vectors. These pBAD33 (Guzman et al. 
1995) and pKK223-3 vectors contain the URA3 selection 
gene and a yeast replication origin obtained from the 
pFL44 vector (Brosius and Holy 1984; Bonneaud et al. 
1991) (see Table 2). Polymerase chain reaction (PCR) 
products were cotransformed with the linearized pBAD33- 
lev or pKK223-3-lev vectors into yeast by Gietz and Schi- 
estl (2007) method. Plasmids were constructed by in vivo 
recombination between PCR products and linearized plas- 
mids in yeast. Candidate vectors were sequenced and 
transformed into the E. coli MC 1061 strain and then into 
the DJ624 Aara. 

Toxicity and antitoxicity assays 

Strains carrying the pBAD33-lev or its derivatives with 
the toxic genes and the pKK223-3-lev or its derivatives 
with the hypothetical antitoxin genes were diluted 10 
times up to a dilution of 10~ 6 CFUs, and 10 fxL was pla- 
ted on minimum media supplied with casamino acids 
0.2%, glycerol 1%, ampicillin, and Cm, and either glucose 
1% or arabinose 1% and IPTG 0.01 mmol/L. 

Results 

Isolation of metagenomic DNA and 
construction of a library 

Metagenomic DNA was extracted from 2 L of filtered 
river water and was sheared before cloning in the fosmid 
vector pCClfos (Epicenter®) as described in Material and 
Methods. Twenty colonies were picked and fosmid DNA 
extracted after a copy number amplification step. Fosmid 
DNAs were submitted to a restriction analysis and each 
clone provided a different pattern, which confirmed that 



the library is representative (results not shown). A total 
of 120,000 clones, each containing an insert of ~35 kb, 
were pooled and screened for antibiotic resistance. A first 
screening of 12,500 clones was done on ampicillin-con- 
taining media (50 /(g/mL), resulting in 20 resistant 
colonies (1/625 clones). DNA was isolated from these 
clones and analyzed by restriction, resulting in two differ- 
ent restriction patterns, suggesting the existence in the 
library of two types of Amp R clones. One representative 
of each clone was selected and plated on media containing 
ceftazidime and piperacillin. One clone showed resistance 
to different classes of jS-lactams (ampicillin, ceftazidime, 
and piperacillin) and shows a high MIC against ampicillin 
with many colonies growing in the inhibition zone, even 
at the highest concentration of antibiotic (Fig. 1). 

Identification of the gene(s) conferring a 
/Mactam resistance 

To determine which gene(s) encode(s) the /^-lactams 
resistance, in vitro transposon mutagenesis of the Amp R 
clone was performed as described in Material and Meth- 
ods. After transformation of E. coli cells, the clones were 
selected on plates containing Cm and kanamycin (to 
select for the fosmid- and transposon-containing bacte- 
ria). After replication on plates containing ampicillin, sev- 
eral Amp-sensitive clones were selected and further 
analyzed. The DNA sequence flanking the transposon was 
obtained using transposon-specific primers as described 
in the Material and Methods section. Three of the inser- 
tions causing the loss of resistance to ampicillin occurred 
in a gene open reading frames (ORF5), which encodes an 
aminotransferase, and the remaining six insertions were 
detected in ORF6, which codes for a /^-lactamase (Fig. 2). 
The transposon insertions in ORF5 could have a polar 



Table 2. Oligonucleotides used in this study, 40 bp regions homologous to the pBAD33lev and pKK223-3lev required for in vivo recombination 
yeast are underlined. 



Primer 



Sequence (5'-3') 



ORF24-for- 
ORF24-rev 
ORF25-for- 
ORF25-rev 
ORF26-for- 
ORF26-rev 
ORF27-for- 
ORF27-rev 
relB93-for- 
relB93-rev- 
relE93-for- 
relE93-rev 



Bad-lev-xba1 
-Bad-lev-pst1 
pkk-lev-ecoR1 
pkk-lev-pst1 
Bad-lev-xba1 
Bad-lev-pst1 
pkk-lev-ecoR1 
pkk-lev-pst1 
lev-ecoR1 
•Iev-pst1 
Iev-xba1 
Iev-pst1 



TTTTTGGGCTAGCGAATTCGAGCTCGGTACCCGGGGATCCTCTAGA AGGAGGgagtaATGAAAAATGAAATAATCGAC 
AAATCTTCTCTCATCCGCCAAAACAGCCAAGCTTGCATGCCTGCAG TTATTTAATCGCATCACTGAG 
TGTGGAATTGTGAGCGGATAACAATTTCACACAGGAAACAGAATTC ATGAAAACTAAAGCACTTTTTG 
CTGAAAATCTTCTCTCATCCGCCAAAACAGCCAAGCTTGGCTGCAGTCATTGTAAATTGTATTCTTC 



TTTTTGGGCTAGCGAATTCGAGCTCGGTACCCGGGGATCCTCTAGA AGGAGGgagtaATGGCAAATTTAGAACCATTC 
AAATCTTCTCTCATCCGCCAAAACAGCCAAGCTTGCATGCCTGCAG TTATCTTCTACCTCTTTGACC 
TGTGGAATTGTGAGCGGATAACAATTTCACACAGGAAACAGAATTC ATGTCAAAATG GATTGATAA 
CTGAAAATCTTCTCTCATCCGCCAAAACAGCCAAGCTTGGCTGCAG CTATTTATC GTTGTTGCTC 
TGTGGAATTGTGAGCGGATAACAATTTCACACAGGAAACAGAATTC ATGAAAACAGAAAACCCCGT 
CTGAAAATCTTCTCTCATCCGCCAAAACAGCCAAGCTTGGCTGCAGTTAGAGCCATTTGGTTAAAAC 



TTTTTGGGCTAGCGAATTCGAGCTCGGTACCCGGGGATCCTCTAGA AGGAGGgagtaATGGCTCTAAATATCCTGTG 
AAATCTTCTCTCATCCGCCAAAACAGCCAAGCTTGCATGCCTGCAGTCAATAATGACTGCGGATATTC 



Primers used to clone into the pBAD33lev vector also include a RBS sequence (AGGAGG) and a spacer sequence (GAGTA) as indicated in bold. 
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Figure 1. Growth of the Amp R clone in the presence of different 
concentrations of ampicillin (Oxoid strips). Although some growth 
inhibition could be observed, there is regrowth of colonies inside the 
inhibition zone. 

effect on the expression of the /^-lactamase gene (ORF6) 
downstream, explaining why they result in ampicillin sen- 
sitivity. Figure 2 shows the gene cluster comprising ORFs 
1-6, whereas all ORFs present in the insert are described 
in Table SI and presented in Figure SI. ORF6 is predicted 
to encode a periplasmic protein containing the 
Pfam 13,354 domain, which corresponds to the serine 
/Mactamase-like superfamily. The closest homolog is a 
/^-lactamase found in Herminiimonas arsenoxydans, a 
/3-proteobacterium (58% identity [169/291]) (Fig. 3). 

A phylogenetic tree shows that the /^-lactamase encoded 
by ORF6 clusters with known class A /^-lactamases 
enzymes, although presenting a different sequence 
(Fig. 4). The ORF6 encoded /^-lactamase clusters with 
other similar class A enzymes present in the genomes of 
/J-proteobacteria. 

Global analysis of the Amp" clone 

To gain information regarding the gene functions associ- 
ated with the detected /^-lactamase, the complete clone 



sequence was obtained (Genbank accession number 
KF033132). After aligning the fragments, a single contig 
of 47,977 nucleotides was obtained. Thirty-eight complete 
ORFs were detected, plus an incomplete ORF39 (Table SI 
and Fig. SI). The first ORF encodes a phosphomannomu- 
tase (phosphoglucomutase), ORF2 encodes a peptidase, 
ORF3 a tryptophanyl-tRNA synthetase, ORF4 a methio- 
nine aminopeptidase, ORF5 an aminotransferase, and 
ORF6 the class A /^-lactamase. ORFs 7-10 follow and are 
transcribed in the same forward orientation and they code 
for hypothetical proteins (Fig. SI). ORFs 11-14 are tran- 
scribed in the reverse orientation and correspond to the 
kdpAFCB genes which encode the components of a potas- 
sium transporting ATPase (Table SI and Fig. SI). The 
ORFs 15 and 16, also transcribed in the reverse orienta- 
tion, do not correspond to known proteins. Interestingly, 
ORF 17, also transcribed in the reverse orientation, corre- 
sponds to the hisC gene encoding a histidinol phosphate 
aminotransferase involved in the biosynthesis of histidine. 
Other histidine biosynthesis genes are also found in the 
contig, namely, ORF21 (hisD) and ORFs 30-34 (hisGBIF), 
also transcribed in the reverse orientation. Between hisC 
and hisD, there are three ORFs (18-20) encoding, respec- 
tively, a RelA/SpoT protein (ORF 18) and two outer mem- 
brane lipoproteins (ORFs 19 and 20). ORF22 encodes a 
ribosome-associated GTPase. ORFs 23-29 code for 
hypothetical proteins or proteins involved in toxin activ- 
ity (ORFs 25 and 28) or in chromosome segregation 
(ORF27). ORF28 product is predicted to be a RelE toxin 
and ORF29 a ribbon-helix-helix protein CopG (see also 
section below). ORFs 34 and 35 are transcribed in the for- 
ward orientation and code for a C4-dicarboxylate 
transporter/tellurite resistance protein and a polyphos- 
phaternucleotide phosphotransferase, respectively. ORFs 
36-39 are transcribed in the reverse orientation and 
code for a cobalamin S-adenosyl methionine transferase 
(ORF36) and different subunits of an ethanolamine trans- 
porter. 

ORFs 26 and 28 encode proteins that inhibit 
E. co// growth 

Two of the ORFs were predicted to encode type II toxins 
(Hayes and Van Melderen 2011). ORF25 was predicted to 



algC Peptidase trpS Met aminopep Aminotransf blaA 

I >r^l >l > l > C73 

ORF i 2 3 4 5 AAA JfcJfcg M. 

Figure 2. Genomic organization of ORFs 1-6: ORF1 (a/gQ encodes a phosphomannomutase, ORF2 a peptidase, ORF3 a tryptophanyl-tRNA 
synthetase, ORF4 a methionine aminopeptidase, ORFB an aminotransferase, and ORF6 a periplasmic /^-lactamase. The black triangles represent 
the transposon insertions resulting in the inactivation of the resistance to ampicillin. 
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1 . VI SSLFKRSLYLLAIMTASANANNLPAVELAKPLPQHSVQKQLAQLESSA 

2 . MPHFSSRRSLLLAAIATPFA-GMAIPAFASDQRISRSAIQQEFSRLEASA 

1 . GGRLGVAAINTANNQRIEYRSEERFPLCSTSKVMAVAAILHNSSTDAKLL 

2 . RGRLGLSAFNTGNGMRVQYRAYERFPVCSTFKTIVAAAILQKSTTDKRLL 



1 . ATNLRYTQQEVAASGYAPIAAKHVATGMTVGAMCQATLNYSDNLAMNELL 

2 . NKRIRYNKDEVEQSGYAPITQKHIADGMSIAELCAATLQYSDNAAANFLM 



1 . KELGGTKAVTAYARKIGDNSFRLDRKEPELNSAIPGDLRDTTTPLAMQNS 

2 . QELGGPAAVTAYARTIGDDTFRLDRWEPELNTAIPGDLRDTSTPAAMGKS 

1 . LQTLVLGNSLAPTQRQQLQEWLKNNTTGSARISAGVPKNWLVGDKTGTGK 

2 . VQRLTLGDALALPQREQLVNWLKGNTTGGKRMQAGVPQDWIVGDKTGTGY 

1 . YGTSNDVGVIWPTGAQPIVLAIYFTQNKPDAKANDQVIADATRIWAEFA 

2 . YGTTNDAGVIWPATGAPIVAAIYFTQNEKEARPRDDI IASATRIWAALG 



2 . KTDAKLKLN 

Figure 3. Alignment of the protein sequences corresponding to 0RF6 
encoded /i-lactamase (1) and Herminiimonas arsenoxydans 
/(-lactamase (2). The identical residues are highlighted in gray and the 
amino acids highlighted in yellow are those found in the class A 
/(-lactamases as described in the results section. 

code for a toxin with a Fic/Doc domain found in Doc 
toxins (Engel et al. 2012) and ORF28 was predicted to 
encode a RelE-like toxin (Overgaard et al. 2009). Both 
toxins are translation inhibitors, although acting at differ- 
ent steps of translation. ORF29, flanking the predicted 
RelE toxin, encodes a predicted RelB antitoxin containing 
a RHH (ribbon-helix-helix) domain. On the other hand, 
the predicted ORF25 Fic/Doc toxin is flanked by the 
small ORF24 and ORF26 with no conserved domains. To 
test whether ORF25 (Fic/Doc) and ORF28 (RelE) encode 
functional toxins, they were cloned and tested in E. coli 
using an inducible expression system as described in 
Material and Methods (Leplae et al. 2011). The potential 
toxic ORFs are cloned in an expression vector under the 
control of an arabinose-inducible promoter, which is 
repressed in the presence of glucose (Guzman et al. 
1995). On the other hand, expression of the predicted 
antitoxins is driven by a pTac promoter, which is induc- 
ible by IPTG. Both plasmids, containing the toxin, and 
the antitoxin gene are compatible. The DJ624Aara E. coli 
strain was transformed by different combinations of con- 
struct and viability was assayed on minimal medium 
plates containing either glucose or arabinose and IPTG to 
induce expression of the different ORFs. Figure 5 shows 
that expression of ORF28 (RelE-like toxin) from the 
pBAD28 construct and the pKK control vector consider- 
ably reduces E. coli viability, whereas expression of 
ORF29 from the pKK29 has no effect on viability. How- 
ever, coexpression of ORF28 and ORF29 restores the 
growth of E. coli. This result confirms that the ORF29- 
ORF28 system is indeed a type II toxin-antitoxin system. 



Expression of the predicted Doc toxin encoded by ORF25 
from the pBAD did not show any toxicity (data not 
show). Interestingly, expression of ORF26 from the pBAD 
vector did inhibit E. coli growth. In an attempt to identify 
a potential antitoxin for the ORF26 toxin, the flanking 
ORFs (ORF25 and ORF27) were cloned in the pKK vec- 
tor and their ability to restore viability was tested. Neither 
of these ORFs was, however, able to counteract ORF26- 
mediated toxicity. ORF24 was also tested and was not 
toxic (data not shown). 

Discussion 

Identification of an active class A |3- 
lactamase 

Our first aim was to detect, by functional screening, the 
expression of antibiotic resistance determinants coming 
from the environment, choosing river water, as it was 
hypothesized that the Zenne river in Brussels, being heavily 
polluted, could be a source of drug resistance genes (Ba- 
quero et al. 2008; Martinez et al. 2009; Allen et al. 2010; 
Garcia- Armisen et al. 2011). The Zenne being as well heav- 
ily polluted by the presence of heavy metals (Garcia-Arm- 
isen et al. 2011), it could also represent a source of both 
metal and antibiotic resistance genes (Baker- Austin et al. 
2006). In a previous study, using the same water samples, 
we obtained 32 cultivable isolates showing resistance to at 
least one antibiotic among the nine tested (Cm, erythro- 
mycin, gentamicin, meropenem, aztreonam, tetracycline, 
amoxicillin, ofloxacin, and ceftazidime) (Garcia- Armisen 
et al. 2011). These cultivable bacteria belonged to different 
phyla (Bacteroidetes and Proteobacteriales) with a large rep- 
resentation of the Pseudomonas genus. However, a 16S 
clone library analysis revealed that there was no good cor- 
relation between the cultivable and the noncultivable bac- 
teria as no Pseudomonas 16S sequence could be detected 
(Garcia- Armisen et al. 2011). This prompted us to look 
for novel resistance determinants by direct functional 
screening of a metagenomic library. The power of func- 
tional metagenomics for the discovery of novel antibiotic 
resistance genes has been previously established by the 
analysis of libraries of soil DNA (Allen et al. 2009; Lang 
et al. 2010; Torres-Cortes et al. 2011) and, more recently, 
from gulls (Martiny et al. 2011). Using this approach, new 
/^-lactamase genes, different from the clinical enzymes, 
could be isolated from an Alaskan soil (Allen et al. 2009) 
and aminoglycoside resistance genes were discovered by 
screening a soil metagenomic library (Riesenfeld et al. 
2004). Although gene expression using E. coli as host is not 
guaranteed (Uchiyama and Miyazaki 2009; Ekkers et al. 
2012), our screening nevertheless resulted in the selection 
of a clone that is able to grow in the presence of high 
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CTXM18 Klebsiella pneumoniae (AF325133) 
CTXM8 Citrobacter amalonaticus (AF1 89721) 
KLUA1 Kluyvera ascorbata (AJ272538) 
KLUC1 Kluyvera cryocrescens (AY026417) 
CTXM1 Escherichia coli (AJ4 16340) 

RAHN1 Rahnella aquatilis (AF338038) 

OXY1 Klebsillia oxytoca (Z30177) 
Sed1 Citrobacter sedlakii (AF321608) 
HugA Proteus penneri (AF324468) 
BestHit3 Photorhabdus asymbiotica (NC012962) 
ORF6 

BestHit4 Variovorax sp. (EJL73433) 
BestHitl Herminiimonas arsenicoxydans (NC009138) 
BestHit2 Janthinobacterium sp. (NC009659) 

NMCA Enterobacter cloacae (Z21956) 

BestHit5 Burkholderia multivorans (NC010086) 
ABPS1 Burkholderia pseudomallei (AF326770) 
AST1 Nocardia asteroides (AF279904) 
A Mycobacterium tuberculosis (U67924) 
ROB1 Haemophilus influenzae (AF0221 14) 



PenP Bacillus subtilis (NC000964) 

BENPEN Bacillus licheniformis (V00093) 

■ L2 Stenotrophomonas maltophilia (AF299368) 



■ SMA1953 Sinorhizobium meliloti (NC003037) 
— GES2 Pseudomonas aeruginosa (AF326355) 



■ HERA2 Escherichia hermannii (AF398334) 
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TEM1 Klebsiella pneumoniae (AF309824) 

PSE4 Pseudomonas aeruginosa (J05162) 

DR0433 Deinococcus radiodurans (NC001263) 

tll21 15b Thermosynechococcus elongatus (NC0041 13) 

CPE 1 184 Clostridium perfringens (NC003366) 

FN1583 Fusobacterium nucleatum (NC003454) 
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■ CFXA2 Prevotella intermedia (AF1 18110) 
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■ CGA1 Chryseobacterium gleum (AF339733) 

— CME2 Chryseobacterium meningosepticum (AF033200) 



► CFB group 



1- 



0.1 substitution per site 

Figure 4. Phylogenetic tree of different class A ^-lactamases (see text for details). 



concentrations of /3-lactam antibiotics. The /^-lactamase 
encoded by ORF6 shows 74% similarity with a class A 
TEM /^-lactamase from the /?-proteobacterium H. arsen- 
oxydans (Muller et al. 2007). The metagenomic /Mactam- 



ase clusters well with A-type enzymes (Fig. 4), suggesting 
that it belongs to this class of /^-lactamase enzymes. Fur- 
thermore, our enzyme also contains conserved residues 
typical of class A enzymes in its active site (Majiduddin 
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Figure 5. ORFs 28 and 29 constitute a RelBE toxin-antitoxin system, 
whereas ORF26 is a lone toxin. Serial dilutions of the DJ624Aara 
containing the pBAD33lev control vector or its derivatives with ORFs 
24, 26, or relE (ORF28) and the control vector pKK223-3lev or its 
derivatives with ORF25, ORF27, or relB (ORF29) were plated on M9 
minimal media containing arabinose (1%) and IPTG (0,1 mmol/L) and 
the appropriate antibiotics. Plates were incubated overnight at 37°C. 

et al. 2002; Bos and Pleiss 2008). These residues, high- 
lighted in Figure 3, are the following: a SXXK motif at 
positions 70-73, a SXN motif at about 130-132, a KT/SG 
motif at positions 234-236, and the f2 loop (Majiduddin 
et al. 2002; Bos and Pleiss 2008). The f2 loop (residues 
161-179) is conserved in all class A /^-lactamases and is 
directly involved in the catalytic reaction of the enzymes 
because it positions the general base Glul66. It also con- 
tains an Arg and an Asp separated by 15 residues, apart 
from each other, which stabilize the loop via the formation 
of a salt bridge, which are also present in our sequence 
(Bos and Pleiss 2008). 

Identification of a relBE toxin-antitoxin 
system 

In this work we demonstrated that ORFs 28 and 29 corre- 
spond to a canonical relBE type II TA system. RelE toxin 
and RelB antitoxin have been described in E. coli, but were 
also identified in Archaea (Li et al. 2008, 2009; Overgaard 
et al. 2009; Shinohara et al. 2010). RelB, the antitoxin, has 
a ribbon-helix-helix structure and it represses the expres- 
sion of the relBE operon using RelE as corepressor (Chris- 
tensen et al. 2001; Christensen and Gerdes 2004; 
Overgaard et al. 2009). Interestingly, ORF18 codes for a 
protein with predicted RelA/SpoT activity. In E. coli and 
other bacteria, amino acid starvation triggers the stringent 
response, which is induced by the production of the alar- 
mone ppGpp, synthesized by the RelA protein (Chatterji 
and Ojha 2001; Sharma and Chatterji 2010). The ppGpp 
molecules associate with the RNA polymerase to inhibit 
the transcription of stable ribosomal RNAs, resulting in 
lowered translational levels (Chatterji and Ojha 2001; 



Sharma and Chatterji 2010). At the same time, the strin- 
gent response also strongly induces the expression of the 
relBE locus, resulting in increased amounts of the mRNA 
translation inhibitor RelE, therefore coupling transcrip- 
tional (RelA) and translational (RelE) inhibition (Chris- 
tensen et al. 2001; Christensen and Gerdes 2004). 
Surprisingly, expression of ORF25 containing a predicted 
FIC/Doc domain did not show any toxicity. The FIC 
domain catalyzes adenylylation and is found among others 
in virulence factors translocated into host cells by patho- 
genic bacteria and in Doc toxins. The consensus motif 
HXFX[(D/E)GNGRXXR was shown to be essential for the 
catalytic activity of Fic proteins (Engel et al. 2012), but 
this motif is degenerated in the putative toxin Doc 
encoded by ORF25, which leads to the hypothesis that this 
degenerated Doc lost its adenylylation activity. Indeed, the 
histidine residue of this motif was shown to be essential 
for Doc toxic activity (Magnuson and Yarmolinsky 1998; 
Garcia-Pino et al. 2008) as well as for FIC proteins (Engel 
et al. 2012). In ORF25, this residue is replaced by a gluta- 
mine (QFFFDGNKRTAR) which may explain why overex- 
pression of ORF25 is not toxic despite its predicted Fic/ 
Doc domain. In contrast, expression of ORF26 was shown 
to be toxic for E. coli. However, this toxin does not appear 
to belong to a classical type II system. Toxicity might be 
due to overproduction and not reflect any functional 
property. 

Histidine biosynthesis genes 

Another striking feature of our sequence contig is the 
presence of different histidine biosynthesis genes, which 
are either clustered (hisFIBG) or separated by other ORFs 
(MsC and hisD). However, the gene coding for HisB, 
which catalyzes the sixth and the eight step of histidine 
biosynthesis pathway, is not present in our contig, mean- 
ing that the probability is high that we did not recover all 
the histidine biosynthesis genes in our fragment. 

From which bacterium is the insert coming 
from? 

The sequenced insert presents several interesting features. 
First, the different ORFs do not show a clear phylogenetic 
affiliation. Eleven of the 39 predicted proteins have best 
hits that correspond to /?-proteobacteria, nine to y-prote- 
obacteria, three to a-proteobacteria, four to <5-proteobac- 
teria, and one to e-proteobacteria. More intriguing is the 
fact that other phyla are represented as well among the 
best hits: ORF10 hypothetical protein best hit is with a 
Calothrix protein (Cyanobacteria), and ORF25 encoded 
protein shows the highest similarity with an Acidobacteria 
protein. The most surprising result of the Blast analysis 
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concerns ORF35 encoding a polyphosphate:nucleotide 
phosphotransferase which shows the highest similarity 
with a protein coming from a methanogenic Archeon 
(Euryarcheota). The histidine biosynthesis genes appear 
also to have homologs in quite different bacteria (a-pro- 
teobacteria for HisC and HisG, y-proteobacteria for HisD 
and HisF, a /J-proteobacterium for HisB, and a bacterium 
belonging to the Bacteroidetes for HisI). 
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Figure SI. Genomic organization of the 39 ORFs present 
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the insert of the Amp R clone. 
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